Thermocapillary valve for droplet production and sorting 
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Droplets are natural candidates for use as microfluidic reactors, if active control of their formation 
and transport can be achieved. We show here that localized heating from a laser can block the motion 
of a water-oil interface, acting as a microfluidic valve for two-phase flows. A theoretical model is 
developed to explain the forces acting on a drop due to thermocapillary flow, predicting a scaling 
law which favors miniaturization. Finally, we show how the laser forcing can be applied to sorting 
drops, thus demonstrating how it may be integrated in complex droplet microfluidic systems. 

PACS numbers: 47.61.-k, 47.55. dm 



Microfluidic droplets have been proposed as microreac- 
tors with the aim to provide high performance tools for 
biochemistry. Individual drops may be viewed as con- 
taining one digital bit of information and the manipula- 
tion of a large number of slightly differing drops would 
allow the testing of a large library of genes rapidly and 
with a small total quantity of material [1]. In microchan- 
nels, drops are produced and transported using a carrier 
fluid and typical channel sizes allow the manipula- 
tion of volumes in the picoliter range. Surfactant in the 
carrier fluid prevents cross-contamination of the drops 
through wall contact or fusion [31, U]. However, while the 
geometry of the microchannel may be used to determine 
the evolution of drops and their contents 0, HI ^ the im- 
plementation of real lab-on-a-chip devices hinges on the 
active control of drop formation and its evolution, which 
remains elusive. 

In this letter, we remedy the situation by demonstrat- 
ing experimentally how a focused laser can provide pre- 
cise control over droplets through the generation of a 
thermocapillary flow. In doing so, we develop the first 
theoretical model of a droplet subjected to localized heat- 
ing, yielding a general understanding of the forces acting 
on the drop and a scaling law which favors miniaturiza- 
tion. A carrier fluid is still used for the formation and 
transport of drops but the effects of geometry are aug- 
mented with a local thermal gradient produced by the 
laser beam, focused through a microscope objective in- 
side the microchannel. 

Indeed, moving drops with heat has been a preoccupa- 
tion of fluid mechanicians since the initial work of Young 
et al [6|. Although originally motivated by micrograv- 
ity conditions where surface effects are dominant [y, , 
microfluidics has opened up a new area where bulk phe- 
nomena are negligible compared to surface effects. Re- 
cently, thermal manipulation of drops or thin films rest- 
ing on a solid substrate has received the attention of the 
microfluidics community either through the embedding 
of electrodes in the solid [8, 9] or through optical tech- 
niques im , [l2| . However, the physical mechanisms 
in the transmission of forces when the liquid touches a 



solid wall are fundamentally different from the case of 
drops suspended in a carrier fluid, away from the bound- 
aries [l3| . The latter case has received little attention de- 
spite the advantages that microchannels offer over open 
geometries. 

Our experimental setup consists of a microchannel fab- 
ricated using soft lithography techniques [l^ . Water and 
oil (Hexadecane + 2% w/w Span 80, a surfactant) are 
pumped into the channel at constant flowrates, Qwater 
and Qoii, using glass syringes and syringe pumps. Chan- 
nel widths are in the range 100 — 500 /im and the height 
h is in the range 25 — 50 /im. Local heating is produced 
by a continuous Argon-Ion laser (wavelength in vacuum 
^Ar+ =514 nm), in the TEMqo mode, focused inside 
the channel through a x5 or xlO microscope objective 
to a beam waist ujq = 5.2 or 2.6 /im, respectively. The 
optical approach can be reconfigured in real-time and it 
allows the manipulation inside small microchannels with 
no special micro-fabrication. The absorption of the laser 
radiation by the aqueous phase is induced by adding 0.1% 
w/w of fluorescein in the water. 

A surprising effect is observed when the water-oil in- 
terface reaches the laser spot: In the cross-shaped mi- 
crochannel of Fig. [H we produce water drops in oil 
through the hydrodynamic focusing technique in which 
two oil flows pinch off water droplets at the intersection 
of the channels. In the absence of the laser, drops of wa- 
ter are produced in a steady fashion and are transported 
along with the oil down the drain channel, as shown in 
Figs. [TJa)-(c). When the laser is illuminated, however, 
the oil-water interface is blocked in place as soon as it 
crosses the beam. While the typical drop pinching time 
is Td ~ 100 ms in the absence of the laser, we find that 
we can block the interface for a time which may be of 
several seconds, as shown in Figs. [Hd)-(f) (see support- 
ing video 1). During the time r^, the drop shedding is 
completely inhibited and the volume in the water tip in- 
creases until the viscous stresses finally break it off. The 
drop thus produced is larger, since it has been "inflated" 
by the water flow. 

In the microchannel shown in Fig. [21 we measured the 
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FIG. 1: Microfluidic valve: In a cross-shaped microchannel, 
the oil flows from the lateral channels and the water enters 
through the central channel, (a)-(c) In the absence of laser 
forcing, drops are shed with a typical break-off time [(b) to 
(c)] of 0.1 s. (d)-(f) When the laser is applied, the interface 
is blocked for several seconds, producing a larger drop. In 
image (e), the evolution of the neck shape is shown through 
a superposition of 100 images (2 s). Exit channel width is 
200 /im. Operating conditions are: Qwater = 0.08 /iL/min, 
Qoii — 0.90 /iL/min, beam power P = 80 mW and beam 
waist (jjQ — 5.2 iim.. 



variation of the blocking time with respect to laser 
power and forcing position. We observe that r5 increases 
approximately linearly with the power, above an ini- 
tial threshold, showing a weak position-dependence of 
the laser spot. Furthermore, the inset of Fig. [2] shows 
that the droplet length L varies linearly with r^, as ex- 
pected from mass conservation at constant water flowrate 
L = Lo-\-ri}Q^3,ter/S, Lq being the droplet length without 
laser and S (125 x 30) /im^ the channel cross section. 
The best linear fit to the data gives an effective water 
flowrate Qwater = 0.028 /iL/min, close to the nominal 
value 0.03 /iL/min, showing that the water flowrate re- 
mains controlled even in presence of the laser forcing. 
Thus, the optical forcing provides a tunable valve which 
provides control over droplet timing and size. Similar 
blocking is observed in a T geometry or if the flows are 
driven at constant pressure. However, the blocking is 
only obtained when the light is absorbed, here by using 
a dye. 

We visualize the convection rolls produced by the heat- 
ing by placing tracer particles in both fluids, as shown in 
Fig. [31(a) for a drop that is blocked in a straight channel. 
For pure liquids, the direction of Marangoni flow along 
the interface is directed from the hot (low surface tension) 
to the cold (high surface tension) regions. However, the 
flows in our experiments point towards the laser along the 
interface, indicating an increase of surface tension with 
temperature. This is consistent with previous studies 
that have shown a linear increase of surface tension with 
temperature in the presence of surfactants [l5|, IH, [l^ • 

One important constraint for practical applications is 
the amplitude of the temperature rise. Since the mate- 
rials used in this study have similar thermal properties 
(thermal diffusivity I^th ~ 10~^ m^s~^, thermal conduc- 
tivity Ath ~ 0.5 Wm~^K~^), we estimate the maximum 
temperature in the flow by modeling the heating pro- 
duced by a laser absorbed in a single fluid phase [18|], 
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FIG. 2: (color online) Dependence of the blocking time rt on 
laser power and position (indicated in the picture) for Qwater = 
0.03 /iL/min and Qoii — 0.1 /iL/min, uq — 2.6 /xm. Inset: 
Rescaled droplet length L/Lq vs. the blocking time (laser 
position •), where Lq is the droplet length without the laser. 
The dashed line is a linear fit, ignoring the outlier at Tb = 
0.75 s. 



assuming thermal diffusion as the only energy transport 
mechanism. Considering the measured optical absorp- 
tion of our water /dye solution ath = 117.9 m~^, and as- 
suming that the temperature 100 /im away is fixed by the 
flowing oil at room temperature, we find AT :^ 12 K for 
the temperature rise at the laser focus for a beam power 
P = 100 mW. The temperature gradient is steep near the 
focus, with the temperature dropping to 5 K at 20 /im 
from the beam spot. However, note that given the typical 
flow velocity {U ^ 1 mm/s) and the characteristic length 
scale over which thermal diffusion occurs (L = 100 /im), 
the thermal Peclet number Pe = UL/Dth is comparable 
to unity. Thus, our calculation overestimates the actual 
overheating. 

The force generated by the convective flow on a droplet 
is investigated through the depth- averaged Stokes equa- 
tions, since our channels have a large width/height aspect 
ratio [l^. The detailed modeling will be discussed in a 
subsequent publication; here we limit ourselves to the 
main features: a circular drop of radius R is considered 
in an infinite domain and the flow due to the Marangoni 
stresses is evaluated. Assuming a parabolic profile in the 
small dimension (h) and introducing a streamfunction for 
the mean velocities in the plane of the channel, the depth 
averaged equations, valid in each fluid, are 



r dr dr dO'^ J \r dr dr dO"^ J ' 

(1) 

where the depth- averaged velocities may be retrieved 



from uq = —dilj/dr and Ur = l/r{dilj/dO). The kine- 
matic boundary conditions at the drop interface (r = R) 
are zero normal velocity and the continuity of the tangen- 
tial velocity. The normal dynamic boundary condition is 
not imposed since the drop is assumed to remain circular, 
which is consistent with our experimental observations, 
Fig. 3a. Finally, the tangential dynamic boundary condi- 
tion, which accounts for the optically-induced Marangoni 
stress, is 



(2) 
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where /ii^2 are the dynamic viscosities and u^''^ are the 
velocities in the drop and the carrier fluid, respectively. 
y = dj/dT is the surface tension to temperature gradi- 
ent, which is positive in our case. 

For simplicity, we approximate the steady state tem- 
perature distribution using a Gaussian form T{x^y) = 
ATexp[-((x - R)'^ + y'^)/'^% where AT is the maxi- 
mum temperature difference between the hot spot and 
the far field and w corresponds to the size of the dif- 
fused hot spot, which is significantly larger than ujq [1^. 
The equations are nondimensionalized using AT as tem- 
perature scale, R as length scale, R^y' /ST as force scale 
and ^^y^^^^ as time scale, the remaining nondimen- 
sional groups being the aspect ratio h/R^ the nondi- 
mensional hot spot size w/R and the viscosity ratio 
Jl = /i2/(/ii + 

A typical predicted flow field solving the above nu- 
merical formulation is shown in Fig. [31(b), in which the 
four recirculation regions are clearly visible. The veloc- 
ity gradients display a separation of scales in the normal 
and tangential directions, as observed from the distance 
between the streamlines in the two directions. Indeed, 
it may be verified that the velocities vary over a typi- 
cal length scale h/R m. the normal direction, while the 
tangential length scale is given by w/R. 

Along with this flow field, we compute the pressure 
field, as well as the normal {aff = 2fl^^) and tangen- 
tial [are = P (^^^ + ^F" ~ ^)) viscous shear stresses 
in the external flow. Their projections on the x axis, 
shown in the inset of Fig. EJc), are then summed and 
integrated along to yield the total dimensionless force 
(F) on the drop. Note that the global x component of the 
force is negative and therefore opposes the transport of 
the drop by the external flow. The y component vanishes 
by symmetry and the integral of the wall friction may be 
shown to be zero since the drop is stationary. Numeri- 
cally computed values of FR/h are shown by the isolated 
points in Fig. [31(c) as a function oiw/R^ for different val- 
ues of the aspect ratio h/R. The points all collapse on a 
single master curve, displaying a nondimensional scaling 
law F (X wh/R^^ for small w/R. 

The dimensional form of the force can be obtained, for 
small h/R and w / R^ by considering the three contribu- 
tions separately and noting that the velocity scale in this 
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FIG. 3: (color online) (a) Overlay of 100 images from a video 
sequence showing the motion of seeding particles near the 
hot spot. Note that the motion along the interface is directed 
towards the hot spot, (b) Streamfunction contours obtained 
from the depth- averaged model described in the text. Dashed 
and continuous contours indicate counterclockwise and clock- 
wise flows, respectively, (c) Rescaled nondimensional force 
FR/h plotted as a function of w/R for various aspect ra- 
tios h/R for fL = 3/4. The straight line corresponds to the 
dimensional scaling derived in the text. The inset shows the 
x-component of the distribution along the azimuthal direction 
of the pressure, normal and tangential shear stresses, where 
the solid circle is the reference zero. Their sum yields the 
total force. Channel width in part (a) is 140 /xm. h/R = 0.2, 
w/R = 0.5 for parts (b) and (c) inset. 



problem is imposed by the Marangoni stress. Using the 
separation of scales along the azimuthal and radial direc- 
tions, Eq. [2] becomes (/ii + /i2) ^ ^ where the 
is understood as an order-of- magnitude scaling. This 
yields the characteristic tangential velocity scale 
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ATj' h 
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(3) 



The force due to the tangential viscous shear is then ob- 
tained by multiplying cFro ^ jii2U/h by sin^ w/R and 
integrating on the portion w x h of the interface. 
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The force due to the normal viscous shear can be shown 
to scale like ~ and is therefore negligible. The 
pressure force, on the other hand, derives from a bal- 
ance between the pressure gradient and the radial second 
derivative of velocity. In the present circular geometry, 
similar scaling arguments yield a law for the contribution 
of the pressure force Fp, which follows the same scaling as 
Ft, resulting in the same scaling law for the total force F. 
A rigorous derivation (to be published elsewhere) yields 
the final form of the force including the prefactor: 

F = -2V^^AT,'^. (5) 

This expression is represented (once non- 
dimensionalized) by the straight line on Fig. ^c) 
and agrees very well with the numerically computed 
values. 

The physical value of the force for a typical experiment 
is estimated by taking /ii = 10~^ Nm~^s (water), fi2 = 
3//1 (hexadecane), and extracting 7' ~ 1 mNm"^K"^ 
from Ref. [16]. This yields a force on the order of 0.1 /iN, 
which is of the same order as the drag force on a drop 
in a large aspect ratio channel [20], thus confirming that 
thermocapillary forcing can indeed account for the block- 
ing. Note that the force we calculate is several orders of 
magnitude larger those generated from electric fields [21| 
or optical tweezers [22]. 

This blocking force may be applied at different loca- 
tions in a microchannel by displacing the laser spot. In 
particular, we demonstrate the sorting of drops, a funda- 
mental operation in the implementation of lab-on-a-chip 
devices. Drops are formed, as above, in a cross-junction 
and arrive at a symmetric bifurcation, carried by the 
continuous phase. In the absence of laser forcing, the 
drops arriving at the bifurcation divide into two equal 
parts [Ej], Fig. [11(a). When the laser is applied, the water- 
oil interface is asymmetrically blocked on the right hand 
side while the left hand side continues to flow, Fig.[4](b). 
Above a critical laser power (approximately 100 mW 
for the present configuration), the drop is blocked long 
enough that it is completely diverted through the left 
hand channel (see supporting video 2). Drops may there- 
fore be sorted by accordingly selecting the laser position. 

In summary, we have experimentally and theoretically 



demonstrated the efficiency of laser-driven blocking of 
water-in-oil drops. The theoretical treatment brings out 
two length scales, h/R and w/R. While h and w can be 
thought of as determining the typical scales for velocity 
variations in the radial and azimuthal directions, R enters 
the force scaling as a local radius of curvature rather than 
the actual size of the drop. It is therefore not surprising 
that the blocking force should increase as R decreases. 
On the other hand, the drag force due to the external flow 
scales as R'^ [20^ , implying that the laser power necessary 
to counterbalance the drag quickly decreases with the 
drop size. This, along with the rapidity of viscous and 
thermal diffusion while thermal inertia is reduced, all lead 





n 


Flow ; 


(b) Laser on 

f 














t 








? 


n 


S s 











FIG. 4: Sorting drops: (a) Without laser forcing, a drop 
at a bifurcation divides into approximately equal daughter 
droplets, (b) When the laser forcing is applied, the drop ad- 
vance in the right hand channel is blocked so the whole drop is 
diverted into the left channel. Main channel width is 200 /xm 
and the operating conditions are Qwater = 0.02 /xL/min, 
Qoii = 0.2 /xL/min, and ujo = 5.2 /xm. 



to laws favorable to miniaturization. The generality of 
the process provides a practical new way for acting on 
individual droplets, at any location, while working inside 
the robust environment of the microchannel. 
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